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0^ I Abstract 

^ ; 

Currently favoured extensions of the Standard Model typically contain 
'flaton fields' defined as fields with large vacuum expectation values (vevs) 

^ ■ and almost flat potentials. If a flaton field is trapped at the origin in the 

. early universe, one expects 'thermal inflation' to take place before it rolls 

^ ■ away to the true vacuum, because the finite-temperature correction to the 

■ potential will hold it at the origin until the temperature falls below 1 TeV or 

so. In the first part of the paper, that expectation is confirmed by an estimate 
of the finite temperature corrections and of the tunneling rate to the true 
vacuum, paying careful attention to the validity of the approximations that 
are used. The second part of the paper considers topological defects which 
may be produced at the end of an era of thermal infiation. If the fiaton fields 
associated with the era are GUT higgs fields, then its end corresponds to 
the GUT phase transition. In that case, the abundance of monopoles and 
of GUT higgs particles will have to be diluted by a second era of thermal 
infiation (or perhaps some non-thermal analogue). Such an era will not 
affect the cosmology of GUT strings, for which the crucial parameter is the 
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string mass per unit length. Because of the flat Higgs potential, the GUT 
symmetry breaking scale required for the strings to be a candidate for the 
origin of large scale structure and the cmb anisotropy is about three times 
bigger than usual, but given the uncertainties it is still compatible with the 
one required by the unification of the Standard Model gauge couphngs. The 
cosmology of textures and of global monopoles is unaffected by the flatness 
of the potential. 



PACS: 98.80.Cq, 98.80.-k, 64.60.-i, 11.27.+d, 12.10.-g, 12.60.Jv, 14.80.Hv 
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1 Introduction 



In presently favoured extensions of the Standard Model, the space of the 
scalar fields contains many directions in which the potential is almost flat 
out to large field values.[] In some of these directions the mass-squared is 
likely to be negative, leading to a nonzero vev whose magnitude will be large 
because the potential is fiat. Fields of this kind, with large vevs and almost 
fiat potentials, have been called fiatons 0, and it has long been recognized 
that they may be cosmologically significant @, ||, ||, ||, 0, ||, ||, [1^, |Tl|, [12[ . [| 
Recently a previously almost M unnoticed aspect of flaton cosmology 



was studied, and termed thermal inflation ||T3|, |14[. Thermal inflation is 
made possible by the flatness of the potential, which near the origin allows 
the (positive) finite temperature contribution to the mass-squared to dom- 
inate the true (negative) mass-squared long after the thermal contribution 
to the energy density has become negligible. If it occurs, thermal inflation 
completely alters the standard cosmology, and one purpose of this paper is 
to examine carefully the features of the flnite temperature potential that 
are supposed to lead to it. Taking into account the validity of the various 
approximations, we shall verify that thermal inflation indeed takes place, 
provided that the fleld is trapped at the origin and that it has sufficiently 
strong interactions with light fields. 

When thermal infiation ends, and the field rolls away from the origin, 
topological defects may be produced just as in the more familiar case of a 
non-fiat potential. The second purpose of this paper is to see what effect the 
flatness of the potential has on these defects, and their cosmology. We are 
particularly interested in the case that the GUT higgs particles are flatons, 
so that the defects produced are monopoles and (depending on the GUT) 

-'^Flat directions occur rather generically in supersymmetric theories as a consequence 
of a non-renormahzation theorem |^ , |^ . The renormaUzabihty requirement impHes that 
there are fiat directions protected to all orders in perturbation theory. The terms allowed 
then have couplings suppressed by some power of the Planck mass. 

^Note the etymology. The term 'flaton' refers to the flat potential, not to inflation. 
Conversely, the familiar word 'inflaton' refers to the field which is slowly rolling during 
ordinary inflation. The term 'flaton' was invented long before it was realized that such a 
field can give rise to a different type of inflation (thermal inflation). 
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cosmic strings, but we shall also study the case of global symmetry breaking. 

Recently it has been pointed out that the interaction of other scalar 
fields with the flaton might cause inflation even without thermalization. In 
that case the situation is more model dependent and we shall not study it in 
this paper, though much of our discussion of topological defects will still be 
applicable. 

Thermal inflation is expected to take place after ordinary (slow-roll) in- 
flation, and to last for fewer e-folds. Such a late epoch of inflation can play a 
welcome role in diluting the abundance of unwanted relics, while leaving unal- 
tered the large-scale density perturbation accounting for the cmb anisotropy 
and large scale structure. (As we shall discuss this remark applies whether 
the perturbation is caused by a vacuum fluctuation during ordinary inflation, 
or by topological defects such as strings or textures.) There can be more than 
one epoch of thermal inflation, so that a second epoch can dilute relics left 
over from a first epoch, and we shall see that this is indeed essential if one is 
to have a GUT phase transition with a flat potential. 

Before getting into more detail we need to be precise about what is meant 
by a 'large' vev, and a potential which is 'almost flat'. In this paper we take 
these terms to be defined with respect to the energy scale 10^ to 10^ GeV, 
which is the scale of supersymmetry breaking as defined by the masses of 
the supersymmetric partners of known particles |]T6| . Thus the vev M of a 
flaton field is much bigger than the above scale, whereas the the energy scale 
I y// 1 1/2 (jggned by the curvature of its potential is taken to be of order that 
scale. In particular the mass m of a flaton particle is taken to be of that 
order. Since the potential as well as its first derivative vanishes at the vev, 
its value Vq at the origin is of order m?M'^. 

Each flaton is either a gauge singlet or a GUT Higgs field. A GUT is 
nowadays considered optional because it is difficult to implement, and can 
perhaps be avoided by appealing to superstring unification near the Planck 
scale. On the other hand the observed Standard Model couplings are com- 
patible with the existence of a GUT, with Higgs vevs of order 10^^ GeV (plus 
possible additional intermediate scale vevs [l^])- If one accepts the existence 
of a GUT then it is attractive to suppose that the Higgs fields are flatons 
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because one might then be able to understand the magnitude of the vev 



without inserting parameters by hand |T^ . 

Whether or not there is a GUT, one expects to find some fiatons which 
are gauge singlets. To understand the possibilities here, note first that each 
flaton will be a complex field.0. If there is a single fiaton field, and it is charged 
under a global f/(l) symmetry, then its vev will spontaneously break that 
symmetry. The prime candidate for a global U{1) symmetry is the Peccei- 
Quinn symmetry associated with the axion, and indeed a fiaton model for 
that symmetry has been proposed With more fiaton fields one can 



have higher global symmetries, as we mention at the end of Section V in 
connection with topological defects. 

However there is no need for a fiaton field to be charged under a continu- 
ous symmetry; on the contrary, it is quite reasonable on the basis of current 
thinking to expect fiaton fields which are charged under at most discrete sym- 
metries. For them, the vacuum manifold consists of discrete points instead 
of an entire circle. The possibilities that we have in mind are of two kinds. 
The first are some or all of the moduh fields [|, |0|, 0, 0, |2|, which 
are ubiquitous in versions of supergravity motivated by the superstring. For 
the present purpose, a modulus may be defined as a field whose potential is 
exactly fiat in the limit of unbroken symmetry. It is widely supposed that 
one or more of the moduli is a fiaton, with a vev of order the reduced Planck 
mass mpi = (SttG)^^/^ = 2.4 x 10^^ GeV.Q (Lower values may be possible 
so that GUT Higgs fields may be moduli ||l8l, but for simplicity, the term 
'moduli' will be used in what follows to denote only the case where the vev 
is of order mpi. It is also possible that moduli have masses much bigger 

■^In a supersymmetric theory all scalar fields are complex since the corresponding left- 
or right-handed spin-half field has two degrees of freedom. 

^The origin of a field is defined to be a fixed point of its symmetries. With all relevant 
fields at the origin the potential has zero slope, where 'relevant' means those fields that 
are coupled to the one in whose direction the slope is being defined. We are assuming that 
this is the situation for the fiatons (or at least that any couplings to other fields displace 
the minimum by an amount small compared with the vev) . In the case of a modulus there 
is more than one fixed point with the separation between fixed points of order mpi. In this 
context the vev may be defined as the distance to the nearest fix point, and the statement 
that the vev is of order mpi simply means that it is not close to any particular fixed point. 
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than 1 TeV, in which case they are not flatons.) Now consider the opposite 
possibihty, of gauge-singlet flatons whose potential is not flat even in the 



limit of supersymmetry Extensions of the Standard Model can eas- 

ily contain such fields, and their potentials will resemble Eq. @ below, with 
couplings perhaps or order 1. As a result, their vev is typically much smaller 
than the GUT scale, perhaps of order 10^ to 10^^ GeV. 

2 Flaton cosmology 

The cosmology of flatons is a rather rapidly moving research area at the mo- 



ment, but we briefly summarize it here following [|14|, |15|, For simplicity 
we generally pretend that there is only one flaton field $, whose potential is 
invariant under a U{1) symmetry so that it depends only on the magnitude 
|$|. We also however point to the differences that can occur in the more re- 
alistic case, where there may be several flaton fields and a higher symmetry 
or else no continuous symmetry at all. 

In the early universe, the effective potential of a flaton field is modified 
by its interaction. If the interaction with light particles is of gravitational 
strength, then the effective potential is expected to receive an additional 
contribution of order ±if^|$p p, |^, If it is somewhat stronger one 

might perhaps expect a contribution ±a^if^|$p with a significantly bigger 
than 1. One expects gravitational strength interactions for moduli with a vev 
of order mpi, whereas one expects interactions stronger by a factor perhaps 
mpi/M for a flaton with smaller vev M. These statements refer to field 
values of order the vev, as opposed to values near the origin which we discuss 
in a moment, and they are based on the idea that the coupling to a light 
particle is suppressed by a factor oc 1/M. 

The weak contributions just described can be of either sign. In contrast, 
if the interaction is strong enough to lead to thermalization then one has the 
finite temperature correction to the potential, which at least near the origin 
gives typically a positive contribution ~ T^|$p to the mass-squared.0 For 

^We will not consider a strong coupling of form SVj ^ — A^$t(j)-^2 ^ji^j-^ \^ ^ which 
would result in a negative temperature-induced mass term SVt '■^ 
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a flaton field $, one expects the thermal contribution to be present in the 
regime |$| < T because the effective mass of particles coupling to the fiaton 
field is of order |$|. Section III below is devoted to a detailed study of the 
thermal contribution. 

As a result of these modifications to the potential, a fiaton in the early 
universe will either be at the origin or will have a value which is large but 
not particularly close to its true vev. We discuss these possibilities in turn. 

2.1 Fiaton field initially at the origin 

Let us suppose that a positive effective mass-squared holds the fiaton field 
at the origin prior to full reheating after ordinary inflation. Then the finite 
temperature mass-squared ~ will typically hold it there until the tem- 
perature falls to a value Tend ~ ~ 10^ to 10^ GeV, where —ml is the 
true mass-squared of the fiaton field at the origin. Only then will it start to 
oscillate about the true vev. On the other hand, the energy density with the 
fiaton field at the origin is 

P = Vo + ^g.T' (1) 

where g^, is the effective number of species.0 Because Tend is so low, the 
first term dominates before the fiaton field rolls away, leading to an era 
of 'thermal' inflation. The era starts when T ~ Vq^'^ ~ (mM)^/^, and it 
ends when T = Tend ~ fn, so the number of e-folds of thermal inflation is 
of order (l/2)ln(M/m) ~ 9 + {l/2)\n{M/10^^GeV). As a result, thermal 
inflation cannot replace ordinary inflation; rather, it takes place if at all after 
ordinary inflation has ended, at the low energy scale Vq^^ ~ (mM)^/^ ~ 
10^ GeV(M/10i° GeVy/\ 

(These type of terms do not destabilize the potential since +A$($'''$)" eventually stabilizes 
it.) 

^This estimate does not apply to non-thermalized particles (in particular moduli) whose 
energy may come to dominate that of the thermalized particles before thermal inflation 
commences. In this case we would have an additional contribution to the energy density 
(the oscillating moduli fields typically behave as massive particles). This would somewhat 
complicate our analysis |l4[ but the qualitative conclusions remain unaffected, so we will 
not consider them. 
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After thermal inflation ends the flaton field starts to oscillate about the 
minimum of the effective potential, so we enter an era of matter domination 
by the flaton particles. (As we discuss in a moment the minimum might 
initially be shifted from the true vacuum, but this effect can be ignored 
because the movement of the minimum will be slow.) If each flaton particle 
decays at the single particle decay rate F, this era ends at the time 
In typical models one estimates F"^ = 1007~^M^/m^ where the numerical 
factor 7 is at most of order 1. At least the bulk of the decay products 
thermalize promptly, so setting F"^ equal to the expansion time where 
H is the Hubble constant we obtain the temperature just after flaton decay 
to be 

where g^^ is the effective number of species at temperature Td, set equal to 
100 for the numerical estimate. 

In fact, the assumption that the flatons corresponding to the oscillation 
decay individually at the one-particle decay rate is not correct, because of 
non-linear relaxation effects of which the most studied example is parametric 
H, ig, |15|, P, H, H (see also 0, |l3). As soon as the flaton 



resonance 

fleld starts to oscillate at the end of thermal inflation, non-linear relaxation 
drains off a signiflcant fraction of the oscillation, or in other words destroys a 
signiflcant fraction of the corresponding flaton particles. They are replaced by 
marginally relativistic scalar particles of whatever species have a signiflcant 
interaction with the flaton fleld, including the flaton fleld itself. 

If nothing happens to the produced scalar particles they will become 
non-relativistic after a few Hubble times, and are expected to decay at their 
one-particle decay rate. If, on the other hand, they thermalize then they 
turn into highly relativistic radiation. At the present time it is not clear 
whether parametric resonance can really create particles which thermalize 
successfully. However, it is clear that the flaton component of the produced 
particles cannot thermalize because here one knows that the interaction is 
too weak. Furthermore, one expects that the energy density of the produced 
flatons will be a signiflcant fraction of the total energy density [E3]. Thus, 
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even if the other produced particles thermahze promptly one expects that a 
significant fraction of non-thermalized energy will remain, and that a signif- 
icant fraction of that energy will be in fiaton particles. 

Any thermalized radiation produced by parametric resonance will red- 
shift away, so independently of the details one expects that a few Hubble 
times after the end of thermal infiation the energy density is dominated 
by non-relativistic scalar particles, including the fiatons and perhaps other 
species. Each species will decay at the single-particle decay rate, so we ex- 
pect eventually to find only the longest-lived species, which dominates the 
energy density until it decays. For simplicity we shall assume in what follows 
that this species is the original fiaton. Then, the upshot of this discussion 
is that despite non-linear effects, the eventual reheat temperature is still the 
temperature calculated from the single-particle decay rate, as given by 
Eq. (I). 

2.2 Fiaton field initially displaced from the origin 

Now suppose that the flaton field has a large value in the early universe. 
In contrast with the first case, the field will not now be in thermal equilib- 
rium because its interaction with light particles will be very weak (because 
otherwise the large value would generate a large mass for the would-be light 
particle). A very crude model for the evolution of the effective potential of a 



flaton fleld $ in the early universe is 0, |2^, |3^ 



^(|$|) = - Mf + aH'^{\^ - %f (3) 

where M is the true vev. The second term represents the effect of interactions 
in the early universe, and for the present purpose we can suppose that $0 ~ 
M is a constant, so that the only time dependence comes from H. 

First suppose that a is of order 1, which is the expected value at least 
for a flaton such as a modulus with gravitational strength interactions (cor- 
responding to M ~ mp;). In that case [jl4|, the minimum of V is practically 



at $0 until the epoch H ^ m, after which it moves quickly to the true vev 
M. The flaton will have settled down to the minimum of V before the end 
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of ordinary inflation, so the conclusion is that if a is of order 1 the flaton 
field starts oscillating about the true vev at the epoch if ~ m, with initial 
amplitude M\ ^ M. 

As soon as the oscillations start, non-linear relaxation effects will convert 
a significant fraction of the energy in the oscillating fiaton field. Provided 
that the oscillation continues to take place about the true vev, then the 
quahtative picture will not be altered, just as we discussed already in the 
case of an oscillation starting out at the origin. However, because the initial 
oscillation energy can now be bigger than the potential energy Vq at the 
origin, there is now the possibility that non-linear relaxation effects restore 



the symmetry, so that the oscillation takes place about the origin [|T5|. In that 
case the fiaton field oscillation might give way to some e-folds of inflation, 
which might be thermal inflation, before the symmetry breaks [jl5[. If it 
occurs, this will happen when the energy density is equal to Vq ~ (mM)^, 
which is much less than the energy density p ~ {mmpiY at the epoch when 
the oscillation starts. 

Finally, consider the case a ^ 1. The evolution of the fiaton field is now 
quite different [^, because the movement of the minimum of the effective 
potential is always slow on the timescale of oscillations around it. As a result 
the fiaton is at all times close to the minimum of its effective potential, and 
the cosmological production of flaton flelds is strongly suppressed. (To be 
precise, analytic BH] and numerical EH] estimates show that the amplitude of 



the oscillation about the minimum is reduced exponentially compared with 
the case a ~ 1.) 



2.3 The case of more than one flaton field 

So far we have considered only a single flaton species. If there are several, 
then some may be initially trapped at the origin and some may be displaced 
from it. Those which are displaced start to oscillate flrst (when the energy 
density is p ~ (mmp;)^, corresponding to a Hubble parameter if ~ m), while 
those trapped at the origin give rise to thermal inflation. If several flaton 
flelds are trapped at the origin, the simplest possibility is that they all roll 
away at more or less the same epoch, which seems reasonable because the 
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negative mass- squared is of the same order for all flatons, and so is the finite- 
temperature correction. If that happens, we essentially recover the case of a 
single flaton. However, the rolling away of the first fiaton (say) might trap 
the remaining flatons at the origin through its interactions. One possibility is 
that the particles produced by non-linear relaxation promptly thermalize, in 
which case one may enter a second era of thermal inflation [|l4l . Alternatively 



the remaining flatons may be trapped by non-linear relaxation ||15[, which 
could also lead to a second epoch of inflation. 

In this paper, we are assuming for simplicity that only genuine flatons can 
roll away from the origin. It could happen that a field with a flat potential, 
and zero true vev, couples to a flaton in such a way that its effective potential 
acquires a nonzero minimum by virtue of the zero value of this flaton. Such 
a field could temporarily act as a flaton, leading to a more complicated 
cosmology p5[ which we shall not consider. 



2.4 The entropy crisis and its solutions 

If a flaton comes to dominate the energy density of the universe, then the 
low 'reheat' temperature Eq. (0) makes it cosmo logically fatal if its vev is 
too large. Indeed, successful nucleosynthesis requires > 10 MeV, corre- 
sponding to M < 10^^ GeV, and thermalization of a stable LSP requires 
Td > 1 GeV corresponding to M < 10^^ GeV. (To have electroweak baryo- 
genesis requires Td > 100 GeV corresponding to M < 10^° GeV, but this 
requirement is not mandatory especially as thermal inflation itself provides 
additional possibilities for baryogenesis [l^, |25[|.) 

This possible problem is a concern for any moduli which are flatons with a 
vev M ~ mp;, and it is also a concern for GUT Higgs flelds if they are flatons 
since the corresponding vev is M ~ 10^^ GeV. It was originally termed the 
'entropy crisis' IQ, and was widely discussed in the 1980's |, ^, |, |10 



O, 12 1 . More recently, focusing on moduli, it has been termed the moduli 



problem [20, pT 



Let us consider the status of the 'moduli' problem, in the light of the 
above snapshot of the current status of flaton cosmology. First, it should be 
emphasised that the problem exists only if the flelds in question are actually 
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flatons, with masses ~ 10^ to 10^ GeV. Assuming that this is the case, the 
status of the problem depends on whether or not the flaton in question is 
initially displaced from the origin. If it is, the problem may not occur at all, 
because the flaton may settle down smoothly to the true vacuum without 
appreciable oscillations; this may well occur for the GUT Higgs though it 



is less likely for the moduli pq , |35[| . If it does occur it can be solved by 
a subsequent bout of thermal inflation, associated with a flaton having a 
cosmologically safe vev M < 10^^ GeV. As we remarked earlier, typical 
extensions of the Standard Model indeed contain such fiatons. 

If, on the other hand, the flaton in question is initially at the origin, then 
thermal inflation will occur and when it ends the flaton will certainly oscillate 
with large amplitude. In that case one needs a second era of thermal inflation 



IJ], again associated with some other flaton which has a cosmologically 
safe vev. At least naively this in turn requires prompt thermalization of a 
significant fraction of the original flaton oscillation energy (with non-linear 
relaxation presumably first converting the energy into marginally relativistic 
particles) and at the present time it is not clear whether this is possible. 
If it is not possible, then we reach the important conclusion that flatons 
with dangerously large vevs, like moduli and GUT Higgs fields, cannot be 
associated with thermal inflation. We return to this issue for GUT Higgs 
fields in Section V. 



3 The 1-loop thermal effective potential 

In this section and the following one, we examine carefully the finite tem- 
perature effects that are supposed to make thermal inflation possible. They 
were first described by Yamamoto in 1986 0, but his account was brief and it 
did not discuss the vahdity of the various approximations employed. To our 
knowledge no subsequent authors have improved on Yamamoto's discussion. 

As already remarked, the potential of the relevant flaton fields might 
be invariant under a continuous symmetry such as a global U{1) or a GUT 
gauge symmetry, or alternatively it might have no continuous symmetry. The 
computations in this paper will be applied to a simple toy model, in which 
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there is a single flaton field with a symmetry. Following [0, |T^, we 
take the potential to be of the form 

r = K,-m^<|.t$ + A^($t$)"^' (4) 

where n is an integer power. Here mo is of order 10^ GeV, and for simplicity 
we have included only the leading non-renormalizable term. The magnitude 
of the coupling is model dependent, but for definiteness one can keep in 
mind a value of order 1. We can parametrize the quantum field $ as follows: 

^= ^ + S(l) + iT], (5) 

where is the classical degree of freedom, which equals to the expectation 
value of a/2($), and S(f) and r] are quantum fluctuating fields with zero ex- 
pectation values. The classical part of the effective potential can be then 
written as 

From this expression we can easily obtain the mass m^ of the flaton particle 
given by ml = (PV/d(p\^ = V2M), the vacuum expectation value M = 
($)o, corresponding to the minimum of V, and the height Vq of the potential 
at the origin which one gets by setting V{M) = 0. The result is 

ml = 2{n+l)ml (7) 

M^(-) ^ j-^^^^f^ml (8) 

Vo = (9) 

One naively expects A„ ~ 1 (n > 1) for the couplings, making M ~ 
10^° GeV. On the other hand, it might be that the Planck scale is replaced 
by the GUT scale if there is one, which effectively increases A„ leading |]14| 



to M ~ 10^. Finally, A„ might be a small mass ratio, which could make 
M S> 10^° so that it becomes the GUT scale (appropriate if is a GUT 
Higgs field) or the Planck scale (appropriate if is a modulus). In this paper 
we focus on the first and second cases. 
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3.1 Coupling to bosons 

Now consider the finite-temperature correction to the effective potential. 
First we will argue that the thermal correction from the $ field itself is 
irrelevant in the sense that it can neither trap the field at the origin nor 
cause a phase transition. Then, in order to cure this we will assume either 
that $ couples to another (real) scalar field via a quartic coupling of form 
or gauge the U(l) by adding the corresponding gauge field. We shall 
also consider the effect of coupling to a spin-half field since we have in mind 
the case of supersymmetry. 

The one-loop thermal correction of a bosonic excitation with the mass m 



to the effective potential has the following generic form ||36[| 

Vlim, T) = ^j^ dxx^ In [l - e-V^^+^^j (10) 

Below we will always quote a one loop result for m. One should keep in mind 
that strictly speaking in order to obtain the one loop thermal correction to 
the effective potential from Eq. ([10|) one should use the tree level value for 
m. In some cases the ring-improved one loop potential (which is a result of 
resumming certain so-called ring diagrams to all orders) is more accurate. It 
corresponds to using the one-loop thermal value for the effective mass only 
for the 'singular' term — m^T/127r in the high temperature expansion. 

The one-loop value for the masses follow: ~ —mQ+gT"^/ 12 is the mass- 
squared of the $-field excitations, and = rn'^(T = + gT'^ / 12+ gcjp' is the 
mass-squared of x- Note that in this approximation does not acquire a 0- 
dependent correction. Indeed, the lowest order correction to proportional 
to (fP' occurs at the {n + l)-loop level and is of order XnT"^ {T / mpiY^ , which 
is minute and hence negligible. We now see from Eq. (p!OD that the finite 
temperature contribution due to the $ field excitations is to an excellent 
approximation a ^-independent shift in the effective potential. This shift is 
nothing but adding a constant to the energy density and hence it should not 
affect the dynamics of the $ field. 

The complete effective potential is now ^(0, T) = V^cj)) + Z^i ^i("^i) ^)) 
where we are summing over all the boson fields. In the regime nii ~ ^/g(t) » 
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T, the finite temperature correction Vi is negligible. (As we remarked earlier, 
one does not in any case expect thermalization in this regime, since particles 
with mass m <^ must couple weakly to (f) [0.) In the opposite regime 



nii ^ ^/g<P ^ the correction is Vi ~ —T'^C^i^j-n^ with (^4 = 7r^/90, which 
is independent of 0. The conclusion is that at large temperatures T ^ M, 
the potential has a single minimum at = 0. (Recall that M is the position 
of the true vacuum, defined by Eq. (|).) As the temperature drops a second 
minimum develops. At a critical temperature defined by ViTf.^ = 0) = 
V^(Tc, = v^M) the two minima are degenerate.^] This temperature is given 
by 

The bump in the effective potential which causes the phase transition is 
solely due to the x field excitations. Even though there are two degenerate 
minima the field will be trapped at = 0, i.e. there will be no phase 
coexistence since the bump is very large and the tunneling rate is minute, 
much smaller than the expansion rate of the Universe. A large bump at 
the critical temperature indicates that the phase transition is strongly first 
order. As the Universe cools down below T^, trapped in the false vacuum 
phase = 0, the false vacuum energy starts dominating. Using Eq. (||), this 
occurs at the temperature T^j ~ Vo7r^/3(y'*C4, at which Thermal Inflation 
begins. The effective number of particle species is equal 3 in our toy 
model, two from the field and one from the x-field. In general (?* will be 
of order 10^ — 10^. Note that, according to the rough estimate Eq. (p, the 
temperature is of the same order of magnitude as the temperature Tti at 
which thermal inflation begins. 

Thermal inflation ends when the Universe supercools sufficiently so that 
the tunneling rate for creation of the true vacuum bubbles becomes compara- 
ble to the expansion rate of the Universe. We will postpone a more detailed 
study of the bubble nucleation to the next section. 



'^Notice that this is not the temperature Tq defined in ||l3|, 14 ; in this paper we are 
denoting that quantity by Tgnd- 



15 



So far we have focussed on the couphng to a spin zero particle x- If the 
flaton is a GUT Higgs field there will also be a coupling to gauge particles. As 
a toy model, consider the Abelian Higgs model, whose Lagrangian density 
is £ = -F^'^F^JA + {D^^)\D^'^) - V, where = + leA^, is 
the vector field, and F^^" is the field strength. The thermal correction to 
the potential has again the form of Eq. ([To|) but with the masses squared 

= — ?tiq + (3e^)T^/12 for both the physical excitation and the unphysical 
Goldstone boson of the 0-field, and m\ = e^cj)^ + (e^/3)T^, m\. = e^cf^ 
for (one) longitudinal and (two) transverse polarized gauge field excitations, 
respectively. At this level some gauge dependence may come into the play, 
but we will not discuss this issue here. The above analysis caries through in 
almost exactly the same way, except that the masses in the thermal potential 
are now uil, ttit and m^. There are 5 degrees of freedom relevant for the 
onset of thermal infiation, so the temperature at which it begins will be 
somewhat lower, as will be the critical temperature. 

What about the validity of the one-loop approximation? The squared 
coupling of the three dimensional theory is = gT for a scalar theory 
and g^ = e^T for a gauged theory. Each loop added costs an extra factor 
gT IM\ recall that ~ ^t,l — e0, so that we expect the one- loop 

approximation to break down when g\ > M, which means 0/T < ^/g for the 
scalar theory and (p/T < e for the gauge theory. In other words, the one-loop 
approximation can be trusted as one approaches the critical temperature up 
to \T — Tc\ ~ gT, e^T . In this case, unlike in the case of the electroweak 
theory, the one-loop treatment is accurate. The reason is that the maximum 
of the effective potential (bump) is located at 

01 ~ ^ (12) 

This should be taken only as an estimate since the bump is located where 
neither the high temperature nor the low temperature expansion is accurate. 
Nevertheless, this estimate tells us something important about the validity of 
the one loop approximation: the bump is located in the region of validity of 
the one-loop approximation (this is, for example, not true in the electroweak 
theory); it is 1/(7 times larger for the scalar theory (1/e^ for the gauge theory) 
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than the value of at which we expect the breakdown of the one-loop ap- 
proximation. We see that, at the local maximum, the 1-loop potential is still 
a reasonable approximation to the true effective potential and this treatment 
of the phase transition is quite accurate. One can also evaluate the potential 
energy density at the bump. To get a feeling for what it is we evaluate it in 
the high temperature limit and find V{(j)i) ~ Vo — T^7r^(l/90 — 1/648). Hence 
we see that the energy density of the bump (T^7r^/648) is much smaller than 
the energy density difference between the false and true vacua (Vq— T^'tt^/QO). 
This suggests that the wall of the critical bubble is thick, which will be used 
in the following section where we study nucleation. 

3.2 Coupling to fermions 

To illustrate what happens in the case of a scalar field coupling with fermions, 
we include in the global U(l) potential a Yukawa coupling of the scalar to a 
spinor field yf(()'4)ip. We will then have a contribution to the effective potential 
of the form 

V^{muT) = -K^ r dxxHn(l + e-V^^^^^) (13) 

271^ Jo \ ) 

where the fermion mass is ^«j(0) = This yields a critical temperature 
very close to the scalar coupling case of Eq. (|TT|) with a minor numerical 
difference: since each fermion contributes to the free energy 7/8 times less 
than a boson, with the same number of degrees of freedom, Tc will be higher 
by a factor (8/7)^''^. The local maximum 0i will also be the same (replacing 
e by ?//). One should note an important difference between this phase tran- 
sition and the electroweak phase transition. Since the main reason for the 
free energy valley at around ~ is not the singular cubic term — mjT/127r 
in the high temperature expansion of Eq. ([T0|) , but the exponential suppres- 
sion of the population density for the mass rrif >> T, the only difference in 
studying the cases of coupling to bosons and fermions is the trivial one of 
counting degrees of freedom. The total effective number of degrees of free- 
dom is g^{tot) = If/* (fermion) -|- (boson). This remark will be also valid 
when studying the nucleation problem in the next section. 
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4 Tunneling 



We are taking the standard route as advocated by Callan, Coleman |3^, and 
Linde in which we assume that the spherical three dimensional action 



-(V0f + n0,T) 



(14) 



determines the nucleation rate per unit time and volume 

Y / C \ „rp2 



. q X 3/2 

V \txTJ 



(15) 



This is an approximate formula but it is sufficiently accurate for our purposes. 



A more accurate treatment of bubble nucleation can be find in |^ . The rate 
in Eq. (|15|) has to be compared with the expansion rate of the Universe 
to determine the rate of bubble formation. More precisely assuming 
a constant bubble expansion velocity, at any given time the fraction of the 
Universe remaining in the symmetric phase is given by 



exp 



■ 47r 

oo O 



(16) 



In order to solve for the nucleation temperature we expand about the value 
when the integrand is about one: S3 = Ss(Tnud) + (t — t'){dS3/dt), where 
d/dt = {dT/dt)d/dT, dT/dt = -TH, and = {^-nG /?,){n^g,{T)T^ /^Q) is 
the Hubble constant. Using a saddle point approximation we find that the 
integral is of order one when 



(HdSs/dT) 



tTTI 



nucl I 



3/2 



vrT, 



nucl 



(17) 



In order to solve this equation we have to find what is the bubble action Ss 
as a function of temperature. 

The solution to the (spherically symmetric, classical) equation of motion 
for 0: 

2d, d 



dr"^ ^ ^ r dr 



-V{<P,T) 



(18) 
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with the boundary conditions: 



^|.=o = 0, 0(r^oo)=O (19) 



specifies the wall profile. The bubble action is 
S's = 47r / dr 



dr 



+ V{<P,T)-V{0,T) 



(20) 



Since Eq. ([18|) is analytically intractable, various approximation schemes have 
been developed. In general when the height of the local maximum is small 
(large) in comparison to the false vacuum energy, the thick (thin) bubble wall 
approximation gives a reasonably accurate answer to S^. First one can show 
using the thin wall approximation, which is accurate around Tc, that the 
three dimensional action: 5*3 ~ (27r/3)o"^/y(0, T)^ is indeed huge; here a = 



j<p,v,v,. ,-v(u,i ^2[V{(j), T) - V{0, T)]d<i) ~ T^/^ is the surface tension of 
the bubble, and near the critical temperature V"(0, T) = ^{T^-T^) « T^. 

When the temperature drops significantly below T^, such that Vbump << 
V{0, T), where Vbump is the difference between the maximum of the potential 
and V{0,T), the thick wall approximation becomes appropriate. Here we 
present a heuristic discussion. The three dimensional action for the critical 
bubble can be approximated as follows: 

Ss ^ 2nR,^' - —Rl[V{0,T) -V{^,T)] (21) 



2^^ {2^J " 

where is the field value at the centre of the bubble. The critical action is 
extremal with respect to R and (p: 

The first condition gives for the size of the critical bubble and the action: 

R^ ^ I S, ^ ^-1 I (24) 

^ 2[V(O,r)-V(0,T)]' 3 J2[V{0,T)-V{<P,T)] ^ ' 
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which can be re-written as 

An ^3 



^9 mo _ 2(^r/mo)2(C4/7r2) + 2(^T/mo)2(y./27r)3/2 exp -y, 

(25) 

The variation of with respect to ip gives then the value of at the origin 
of the bubble: 



^^-^^) (26) 

ttM mo / 



2 



which is obtained by neglecting the exponentially small term in Eq. (psf ) so 
that the critical action is 



T TT Vmo/ 



We can now use this result and Eq. (|T6|) - ([T^) to obtain the nucleation 
temperature Tnud- Q 

Snuci = — — — 6 + 41n — — - + - In — ~ 17 + 4 In — 28) 

where we assumed that the bubble is supersonic: f ~ 1; this is plausible 
since supercooling is very large, temperature is about Tnud — 2mQ. For 
somewhat different values of parameters Tnud changes only slightly. To a very 
good approximation Tnud is proportional to mo- (Here we have neglected a 
logarithmic dependence of the form (Inmo)^^^.) 

These results are in substantial agreement with that of Yamamoto [0], 
who found that the bubble action is (up to In mo) proportional to T'^/rriQ. 
What we have done is to look more carefully at the somewhat delicate as- 
sumptions which have to be made in order to derive it. We have found that 
at a reasonable level of confidence they can indeed be justified by carefully 
considering the regime of validity of the various approximations. 

*Note that Tnud for all practical purposes coincides with the temperature at the end 
of thermal inflation Tend, discussed in section 2 above. 
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5 Topological defects with a flat potential 



Now we come to the second part of our investigation, which is to look at 
topological defects associated with a flat potential. In the case of a single 
(complex) fiaton field with a U{1) symmetry, which we have focussed on 
for the sake of simplicity, these are cosmic strings. If instead there is only 
a discrete Z„ symmetry they are domain walls. Finally, if there are two 
or more fiaton fields there could be a higher continuous symmetry, which, 
depending on the pattern of symmetry breaking, could lead to monopoles 
and/or textures (or to no defects at all). We shall consider all these cases. 

In order for the defects to form, the symmetry has to be restored at early 
times, which as discussed in Section 2 may or may not happen for a given 
fiaton field. Assuming that it is indeed restored the defects form when the 
field rolls away from the origin, and since we are considering a fiaton field 
this occurs at the end of an era of thermal inflation. Thus, we are considering 
the formation of topological defects at the end of such an era. 

We shall focus particularly on the case where there is a GUT whose Higgs 
fields have a flat potential, schematically of the form Eq. (^. As mentioned 
already such a flat potential seems quite natural, because one can then hope 
to generate the required vev from scales already present in the theory [|T^]. 
For example, in Eq. (^, the coupling A„ might be of the form {nig/mpiy, 
where rrig is related to vevs arising from gaugino condensation, and p is a 
positive integer. This paradigm is of course very different from the usual one, 
where the potential of the GUT Higgs is supposed to be non-flat so that it 
has the form 

V = X{<^^<^ - M'^ f (29) 

with A ~ 1. Such a form is indeed natural for the Standard Model Higgs 
where M is only of order the susy breaking scale 10^ to 10'^ GeV, but it seems 
far less natural for a GUT. In other words, there seems to be no reason to 
suppose that the vev of a GUT Higgs is of the same order of magnitude as 
its mass. 

GUT symmetry breaking with a flat potential is completely different from 
the usually considered case. The GUT phase-transition occurs at the end of 
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an era of thermal inflation. Just before it occurs, the potential energy density 
Vo of the Higgs field accounts for most of the energy density in the universe, 
and after the transition this energy is initially converted into a homogeneous 
oscillation of the Higgs field (away from any topological defects), correspond- 
ing to extremely non-relativistic Higgs particles. Through non-linear relax- 
ation effects this initial era is quickly followed by one in which at least a 
significant fraction of the energy density resides in marginally relativistic 
scalar particles, of various species including the Higgs particles. There may 
also be a significant fraction in thermalized (and hence highly relativistic) 
particles, but this question has not been settled at the present time. 

Because the potential is fiat, the GUT Higgs particles are light (mass m ~ 
10^ to lO^GeV), and according to Eq. (H) they decay after nucleosynthesis. 
In order not to upset nucleosynthesis their abundance must be diluted by 
a second era of thermal inflation.^ Such an era requires the thermalization 
of a substantial fraction of the Higgs particles, and, taking the fraction to 
be of order 1, but not very close to it, the abundance of the Higgs particles 
immediately after thermalization is n/s ~ Vq^^ / {gi^^m) ~ 10^. From pi|] , 
their abundance at the epoch of nucleosynthesis must satisfy n/s < 10~^^, so 
the second bout of thermal inflation should dilute the abundance by a factor 
of at least 10^^ 



The actual dilution factor actually provided by thermal inflation is [|T^| 



In this expression m and M refer to the second era of thermal inflation, as 
do Tend (the temperature at the end of thermal inflation) and To (the reheat 
temperature). The quantity e < 1 is the fraction of energy left unthermalized 
after non-linear relaxation effects have occurred after the end of the second 
bout of thermal inflation, which we set equal to 1. 

In making the above estimates we have set Vq = m^M^ and Tend = rn = 
10^'^ GeV for both eras of thermal inflation. The expected range for the vev 
of the inflaton causing the second bout of thermal inflation is 10^ GeV < M < 



^Or perhaps a non-thermal analogue |E5[, but we shall not consider that possibility. 
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10^^ GeV (with the lower hmit coming from the discussion after Eq. and 
the upper hmit from the requirement of decay before nucleosynthesis) . Since 
several parameters have been set equal to fiducial values, a second bout of 
thermal inflation seems well able to provide the necessary dilution. The only 
problem is to ensure that it occurs, bearing in mind the fact mentioned earlier 
that the possibility of thermalization immediately after thermal inflation ends 
has not been demonstrated. 

Assuming that the GUT particles can be somehow diluted, we proceed 
to a discussion of the defects produced, starting with monopoles. 

5.1 Monopoles 

It is well known that monopoles form when the vacuum manifold M of a 
gauge group G contains non-shrinkable 2 surfaces, i. e. when 7r2(M) 7^ /, 
where 112 is the second homotopy group. In the light of a theorem in homo- 
topy theory, this may be rephrased as: monopoles form whenever a grand 
unified semisimple group G breaks down to a group H which has nontrivial 
fundamental homotopy group vri(II) 7^ J, i.e. H contains non-contractable 
loops. A simple example is an H which contains at least one U(l) as a factor, 
just like the Standard Model! It is argued in |^ that since supersymmetry 
does not change the connected structure of the (super)-Lie group, the same 
consideration as above is valid for a supersymmetric GUT. Since rather gener- 
ically the symmetry of a low energy theory contains non-contractable loops, 
one must check what would be the cosmological implications of monopoles' 
formation, i.e. whether there is a monopole problem. 

It is convenient to specify the monopole abundance rzmon as a fraction of 
the entropy density s, because rimon/s is constant as long as monopoles are 
not destroyed, except when particle decay (or some other non-equilibrium 
process) increases the entropy per comoving volume. Consider first the ob- 
servational upper bound on the present monopole abundance. The mass of 
a gauge monopole for a given vev M is not much affected by the flatness 
of the potential, because it is known that the mass goes to a finite limit 
as the potential becomes flatter and flatter (Bogomolny bound) This 
leads to a simple and reliable bound on the present monopole abundance of 
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nmon/ s < 10^^^, coming from the requirement that the monopole density 
be less than criticaL This bound is stronger than that derived from the re- 
quirement that monopoles do not affect nucleosynthesis. The next bound 
one can derive is based on the requirement that the galactic magnetic field 
not be dissipated through acceleration of monopoles from the time of galaxy 
formation |Q. Strictly speaking this bound applies only to the monopole 
abundance in galaxies. The original Parker bound is rimon/s < 10~^^. This 
was argued to be too strong; when the depletion of monopoles, as they get 
accelerated and ejected from the galaxy, is taken into account [^, the bound 
becomes Umon/s < 10~^°. A stronger bound was found using the fact that 
a grand unified monopole catalyzes nucleon decay when captured by a 
neutron star; the bound is Umon/s < 10~^^ A more stringent bound 

can be derived based on monopole capture and consequent catalytic nucleon 
decays in the main sequence stars Umon/s < 10~^^. The flatness of the 



potential is not expected to alter this bound because the strong interaction 
cross section for the monopole is independent of its radius ^ . 

In order to find out whether there is a monopole problem, we ought to 
compare the above bounds with the amount of monopoles that form at the 
phase transition. Since we need a second epoch of thermal inflation (to 
dilute the GUT Higgs particles), a significant fraction of the energy density 
had better thermalize immediately after the GUT Higgs phase transition, 
leading to a temperature Tqut ~ (mM)^/^ ~ 10^ GeV. (Here m ~ 10^ to 
10^ GeV is the mass of a typical GUT Higgs particle with a fiat potential, and 
M ~ 10^^ GeV is the vev of a typical GUT Higgs field.) This is much bigger 
than the temperature T ~ m just before thermalization, but it is still far 
below the temperature T ~ M which one would obtain after thermalization 
with a non-flat potential. In the case of a non-flat potential the temperature is 
high enough to initially annihilate the monopoles, so that when annihilation 
ceases there is of order one monopole per horizon volume irrespective of the 
initial abundance. For a flat potential this is not the case, and we have to 
consider the initial abundance. 

A lower limit on it is provided by the causality bound which states 
that on average about p ^ 1/8 monopoles form per causality volume Vc — 
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(47r/3)iJ^^, where H^^ is the causahty radius. The probabihty p ~ 1/8 can 
be obtained from an argument based on the Kibble mechanism: the field is 
uncorrelated on superhorizon scales, and the probability that the randomly 
oriented field in different correlation volumes will cover more than half of the 
vacuum manifold is about p = 1/8. This bound gives an initial monopole 
density rimon — {S/4:Tt)pH^, corresponding to 



nmon\ . 8p (moM)3/2 27 /Q1^ 

^ 1/4 ^3 10 (31) 

g* ^pi 



We have used for the entropy density s ~ gy^el^^^{'mQM)^^'^; here M denotes 
the GUT Higgs vev, mo the curvature of the potential at the origin and 
^conv is the fraction of the energy in the field that decays and thermalizes. 
The numerical estimate in ([31| ) is obtained by setting rriQ ~ lO'^GeV, and 
M ~ lO^'^GeV. 

To obtain an actual estimate, remember that we have argued that the Uni- 
verse undergoes a strong supercooling and the phase transition is strongly 
first order proceeding via bubble nucleation. An estimate of the typical bub- 
ble size at collisions will allow us to estimate monopole production. Since 
each bubble nucleates with a random phase at the vacuum manifold, when 
four bubbles collide a monopole (or an anti-monopole) forms with the prob- 
ability p ~ 1/8. The only missing information is now the typical bub- 
ble size when they collide. A typical bubble nucleates at the temperature 
Tnud — 2mo, as specified by (^TD and Eq. (^8]). The typical bubble size can be 
obtained from Eq. (plGf) as follows. The integral in the exponent is dominated 
by bubbles of size 

/ o 1 3 T 1 11 

Hence, the typical bubble, when it collides, is about l/20th of the horizon 
size 1/H. This means that we have an estimate for monopole production 
which is about 20^ ~ 10^ times above the causality bound in Eq. (|3TD: 

nmon\ _ ^Q_23 (33) 

bubble 



25 



According to this estimate the dilution factor provided by the second bout 
of thermal inflation is ample. 

In the discussion so far we have supposed that the monopoles are isolated, 
but in many cases they are connected by cosmic strings |^9|, If this is 
the case then strong effective forces will be induced by connecting monopoles 
and anti-monopoles with strings which will cause efficient annihilation all 
the way to the horizon scales, establishing a scaling solution. (In this respect 
the cosmology of the strings connecting monopoles resembles the cosmology 
of global monopoles discussed below. This does not mean that monopoles 
connected with strings may affect structure formation, as it is in the case 
of global monopoles discussed below, since the ordering field dynamics of 
global monopoles is very different than that of the magnetic monopoles. In 
these models, if any cosmologically interesting effects exist, they would be 
produced by the strings.) On scales larger than the horizon causality prevents 
annihilation. Nevertheless, this annihilation mechanism will suffice to keep 
the monopole density at the level of a few per horizon at all times and hence 



solve the monopole problem EOl, H3 



On the assumption that there is no monopole problem, let us proceed to 
a discussion of cosmic strings. 



5.2 Cosmic strings 

Cosmic strings form in the case when the vacuum manifold M of the broken 
theory contains non-shrinkable loops so that 7ri(M) 7^ /. In other words 
(using some homotopy theory), they form if a semisimple grand unified group 
G breaks down to a group H which is not connected, i.e. 7ro(H) 7^ I. (Recall 
that connectedness of a group is not changed by supersymmetry.) Since the 
symmetry groups of low energy theories are often connected, cosmic strings 
are a less generic feature than monopoles, but they do occur for some choices 
of the GUT gauge group ||4^. (As an example, SO(IO) may break down to 
the Standard Model enhanced by an additional matter parity Z2, which is, 
of course, not connected.) 

On scales below the horizon, the cosmic string network has a 'scaling' 
configuration which is more or less independent of its previous history. As a 
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result, the cosmology of the strings will not be significantly affected by the 
ten or so e-folds of inflation taking place after their formation. Instead, it 
will depend mainly on the string mass per unit length /i. As is well known, 
the strings are candidates for the origin of both the cmb anisotropy and large 
scale structure, with some value /i ~ (10^^ GeV)^. If this turns out to be the 
case then /i will be determined, and in any event one already has an upper 
bound on fi. In contrast with the more commonly considered model of a 
smooth, Gaussian primordial perturbation, the cosmic string prediction has 
yet to be accurately calculated over a wide range of scales, but as a guide we 
may take the prediction for the low multipoles of the cmb anisotropy. Here 
it is known ^ 0] /i ~ 2 x (10^^ GeV)^, with an uncertainty of perhaps 
a factor 2 either way.|3 If the effect is observed and is due solely to strings, 
one will therefore deduce a value of fi somewhere in the range 

/i ^ 1 to 4 X (lO^^GeV)^ (34) 

which imposes an upper limit 

/i < 4 X (lO^'^GeV)^ (35) 

What do these results for /i tell us about the scale of GUT symmetry 
breaking, as defined by the vev M? With a non-fiat potential of the form 
Eq. (^) it is known |5l|] that fi ~ vrM^. Using this result, an observed signal 
in the cmb anisotropy will imply a value 

5 X 10^^ GeV < Af < 1 X 10^^ GeV (36) 

and one already has an upper limit 

M < 1 X 10^^ GeV (37) 

We want to know how much this estimate is changed by the flatness of the 
potential, and in contrast with the case of the monopole mass the question is 

^°One of the uncertainties in this estimate is the fact that smaU-scale structure in the 
string network breaks the equahty of the mass-per-unit length and the tension, by a few 
tens of percent [ p3[ . The estimates being referred to, on the other hand, are performed 
under the assumption that these quantities have a common value fi. 
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non-trivial because /i does not tend to a finite limit as the potential becomes 
flatter and flatter. 

We work with the same toy model that we considered earlier, namely the 
case of a U(l) gauge symmetry and a potential of the form (|]) with n = 2 (0® 
coupling). The Lagrangian density of the flaton coupled to the U(l) gauge 
field is 

C = -^F^uF^" + |D^<I>|2 - V{^) (38) 

where = + ieA^j,, F^^ = d^Ay — duA^ and V is the tree-level, zero 
temperature scalar potential. The equations of motion in the Lorentz gauge 
{d^A^ = 0) for this Lagrangian are 

+ ^ = (39) 
d^F^'"' + ie{^*D^^ - (D'^$)*$) = (40) 

We now make a (static) cylindrical ansatz for $ and A^ (setting, for simplic- 
ity, the winding number to one). 

^ = M f{mvp)e'^ and A=—a{mvp)0 (41) 

ep 

where (p, 9, z) are the cylindrical coordinates, / and a generic functions, M 
is the vev of the $ field (see equation ||), and my = eM is the vector particle 
mass. The equations of motion will then read: 

r + - ^/(« - 1)' - Pfif - 1) = (42) 
a" - ^ + 2/2(1 - a) = (43) 

where we have defined ^ = myp and /? = ml /my. 

The asymptotic behaviour of the solutions to these equations are similar 
to the ones for the 0^ potential 

f M, f ao^^ as ^ ^ 0; 

^^11- M-'^'exp{-y/^0, 1 1 - a,e^'exp{-0, as ^ ^ oo. 
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We now are interested estimating the energy per unit length, /i = / edV, 
of these strings. The integral is over all space, and e is the energy density, 
which in this case reads as 

2 



9$ 


2 


dp 


+ 



19$ 
'p'de 



ieAg^ 



(44) 



where i? = V A A is the magnetic field. With our ansatz, the energy per unit 
length becomes 







(45) 



In the usual case, with /3 ~ 1 and the quartic coupling, the integral is close 
to one (it can actually be shown that it is one when (3 = 1) and the mass 



scale of the string is given by the vacuum vev [^, ^ . One expects a similar 
result here, but our coupling is much smaller than usual with (3 ~ 10"^''^. It 
is not very difficult to estimate the energy per unit length and the profile of 
these static vortices numerically. Figure 1 shows the results thus obtained 
for a large range of values of the (3 parameter for both the quartic and cf)^ 
cases. The energy decreases logarithmically with small (3 and so the very 
small coupling changes the value of p by only a small factor of two orders of 
magnitude. The difference between the two types of potential is also seen to 
be small (around 20% for /3 ~ 1) and getting smaller with decreasing (3. 

that when (3 is very small, p becomes rather insensitive to the value of /3, 
with a value p ~ (0.09 — 0.14) x (to be compared with the value p ~ vrM^ 
for a non-flat potential with (3^1). This estimate applies in a rather large 
range of the coupling (3 ~ 10^^'^ — 10^'^^. It follows that an observed signal 
in the cmb, corresponding to the range deflned by Eq. (0), will correspond 
to 

3 X 10^^ GeV < M < 6.5 X 10^^ GeV (46) 
which converts to an upper bound on M of 



M < 6.5 X 10^^ GeV 



(47) 



For a given GUT theory, the scale Mqut at which the three Standard 
Model couplings unify is determined by their values measured at low energies. 
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Beta 

Figure 1: ^ as a function of j3 
The energy per unit length of the static vortex solution for the cp"^ potential (full 
line) and the (j)^ potential (dashed line). The energy per unit length is in units of 

In simple GUT models, Mgut is found to be around 2 x 10^^ GeV, and 
with a single Higgs field one would have gM ~ Mqut where g is the gauge 
coupling at the unification point. This coupling in turn is typically given by 
g'^/{4:7i) ~ 1/25, corresponding to 5^ ~ 1, so that with a single Higgs field 
one would have M ~ Mqut- However, in a realistic model there will be 
more than one GUT Higgs field and all of the analysis in this subsection will 
require modification. 

Taking into account the present uncertainties, the GUT unification scale 
required to make strings a candidate for the origin of large scale structure is 
consistent with the one deduced from low-energy data, both for a fiat and for 
a non-fiat potential. In particular, the upper limit on the unification scale 
deduced from the cmb anisotropy be no bigger than that observed by COBE 
is consistent for both forms of the potential. On the other hand, it is clear 
that in the future one or other form of the potential might be preferred, since 
we have seen (Fig. 1) that the ratio /x/M^ is significantly different for the 
two cases. 



30 



5.3 Global symmetry defects 



Although a GUT is initially constructed to satisfy a gauge symmetry, global 
symmetries could occur accidentally and be spontaneously broken by fields 
with vevs at the same scale ~ 10^^ GeV. Global symmetries, notably the 
Peccei-Quinn U{1) symmetry have also been proposed at lower energies. In 
addition to the possibilities of monopoles and strings, there are now the 
possibilities of domain walls and textures. We discuss textures first. 



5.3.1 Textures 

It is known that textures with a vev M ~ lO"*^^ GeV can significantly affect the 
cmb anisotropy and large scale structure [Q. Again, the central question we 
want to answer is whether the flatness of the potential alters this situation. 
As with strings, the configuration approaches a scaling solution so the only 
thing we have to worry about is the energy of the textures for a given vev 
M. The scaling solution for textures is characterized by continuous texture 
formation on the horizon scales and their subsequent collapse, as governed 
by the scalar field ordering dynamics. 

It turns out that the relevant mass scale is not affected by the flatness 
of the potential. We will now argue that this is indeed so. The size of the 
self-coupling A does not matter as long as it is large enough to conflne $ 
to the vacuum given by |<I>| = M. Since we are interested in the epoch 
relevant for structure formation, the relevant scale is T^ec ~ leV, which is 
much below any scale in the microscopic theory. Indeed the 'flatness scale' 
is given by the curvature of the potential at its minimum m = V"{MY^'^ ~ 
mo ~ 1 TeV >> Tdec- A more rigorous argument can be constructed by 
looking at the evolution equation for textures in an expanding Universe. 
This equation has an effective self-coupling re-scaled as Xejf ~ (a(T)/ao)^A 
\^% , where a is the scale factor of the Universe. At T ~ Tdec, the ratio 
a{Tdec)/ao ~ Tcur/Tdec is so big that Ae// ~ {rrio/Tdec)'^ » 1 is still very 
large. This means that the effective coupling for all practical purposes is 
large and the fleld is conflned to I'I'I = M. This implies that the non-linear 



cr-model is a very accurate description of the texture dynamics, even for 
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flat potentials. Normalising to the 10° scales reported by COBE, and using 
the cosmo logical parameters = 1, A = and h = 0.5 the result may be 
read off from M ~ 1.4 x lO^^GeV. 

A possible realization of a texture model within grand unified theories 
is proposed in in which the GUT of the form G x SU{3)fam (where 



SU{3)fam is an additional global family symmetry) has been considered. A 
supersymmetric GUT group with a flat Higgs potential may have a similar 
form. The family symmetry SU (3) jam breaks with the GUT, so it is plausible 
to assume that the texture potential is also flat. 

5.3.2 Global monopoles and domain walls 

The dynamics of global field ordering with SO (3) symmetry (global monopoles) 
can be studied using similar numerical techniques as developed for textures 
[0 . Also the ordering dynamics can be to a good approximation represented 
by a non-linear cr-model , [Q. They eventually reach a scaling solution 



in which 'monopoles' form on the scales of the horizon and then collapse to 
the core size given by ~ I/ttiq. Since the interactions between monopoles 
are strong and long range, there will be a strong tendency for monopole- 
anti-monopole annihilation which will keep the number of monopoles (anti- 
monopoles) at the level of order one per horizon at any time. In that re- 
spect, the monopole dynamics resembles the dynamics of textures, with an 
important difference of an occasional monopole - anti-monopole annihila- 
tion event, which will tend to imprint a strong non- Gaussian signal on the 
cmbr. Recalling the assumptions made in their calculation concerning the 



values of the cosmological parameters we have read off from |£2[ that global 
monopoles have an observable signature on large scale structure and/or the 
cmb anisotropy for some vev of order M ~ 1 x lO^^GeV. 

For walls we can make a simple estimate based on the energy per unit area 
of the wall: cr ~ / dz[{d(j) / dzY / 2 + V{z)]. Making the simple substitutions 
d/dz ~ 1//, /rf^ ~ /, V{z) ~ I/q ~ (moM)2, we have a ~ M^/l + VqI, which 
is extremized when / ~ M/y/Vo ~ ^/ttlq so that the energy per unit area 
is about a ~ M^{mo/M), where rriQ ~ 10^ GeV. Compared with the usual 
case one sees that walls are much fatter, but less massive. However, the 
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reduction in mass is far too small to allow them to be cosmologically viable. 
Instead they had better either not form or disappear harmlessly, which is 
quite possible since the discrete symmetry leading to their appearance will 
probably be slightly broken if it is not a gauge symmetry [O]. 



6 Conclusion 

The topics discussed in this paper relate to a considerable body of ongoing 
research at the moment, whose common theme is the form of the effective 
potential in the early universe, and the consequent cosmological effects of the 
scalar fields. Some of this research continues to address the traditional, and 
still very important, question of how to implement an early era of ordinary 
(slow-roll) inflation. Here we have focussed instead on aspects of 'thermal 
inflation', which occurs if at all after ordinary inflation and lasts for only a 
few e-folds. 

The fields which can give rise to thermal inflation are 'flaton' fields. By 
definition, these fields have 'flat' potentials and 'large' vevs, where these 
terms refer to the mass scale 10^ to 10^ GeV. They arise naturally in currently 
favoured extensions of the Standard Model, along with the opposite case of 
fields with 'flat' potentials but zero vevs. 

A flaton field gives rise to thermal inflation if it is trapped at the origin 
by virtue of the finite temperature correction to its effective potential, after 
the thermal contribution to the energy density has become comparatively 
small. According to a calculation made a decade ago by Yamamoto, this 
trapping is expected to occur for a flat potential (provided that the flaton 
field is indeed in the vicinity of the origin), with thermal inflation ending 
only when the local minimum of the effective potential at the origin almost 
completely disappears. Part of our objective was to investigate carefully 
the rather delicate assumptions needed to arrive at this conclusion, using 
some modern perspectives particularly on the tunneling rate out of the false 
vacuum. Happily, we confirm the conclusion. 

The actual cosmology of fiaton fields, which determines whether a given 
fiaton field will actually find itself in the vicinity of the origin so that thermal 
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inflation can take place, is at present a rather rapidly moving research area. 
We have not attempted any new advance here, contenting ourselves with 
a brief snapshot of the situation as it stands at present. Instead we have 
addressed a question which has received relatively little attention, which 
is the nature and cosmology of topological defects forming at the end of 
thermal inflation. One would expect a priori that their cosmology might 
be significantly affected by the flatness of the potential. The answer to this 
question turns out to depend on the type of defect. 

A particular focus for the second part of our investigation has been the 
possibility that there is a GUT, whose Higgs fields have a fiat potential. This 
possibility is quite natural, and it leads to a cosmology very different from 
the usually considered case of a non-flat Higgs potential. GUT symmetry 
breaking, if it occurs, will be preceded by an era of thermal inflation. The 
Higgs particles produced at the transition are cosmologically dangerous be- 
cause they are hght and long-hved, and to dilute them one needs a second 
era of thermal inflation. This period of thermal inflation will also dilute 
the monopole abundance sufficiently, but will not eliminate cosmic strings 
because it lasts only a few e-folds. A simple picture of a string network evo- 
lution in an inflationary Universe is as follows. First as the Universe expands 
the strings will quickly reach the density of about one (long) string per hori- 
zon volume. After that the network freezes out since there exists no causal 
process that could incite nontrivial dynamics. This means that from then 
on the average correlation length grows exponentially as the expansion rate, 
i.e. BS ^{t) — ^{to) e^^^~^°\ After about 10 e- foldings the correlation length 
has not grown more than ^ ~ e^^/H; this scale will come back within the 
horizon after ln(a/ao) ~ 10 which in radiation era means at the temperature 
T = Td X e~^° ~ To X 10~^'^; after this strings will quickly reach the radi- 
ation era scaling solution. This temperature is far above the relevant scale 
for onset of cmbr anisotropics and structure formation, keeping the strings 
cosmologically interesting. 

With all this in mind, we expect that cosmological relevance of cosmic 
strings is specified solely by the mass per unit length of strings. We find that 
the GUT Higgs vev needed for the strings to give an observable signature 
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in the cmb anisotropy is a few times bigger than in the case of a non-flat 
potential, but still compatible with the typical estimate ~ 2 x 10^^ GeV of the 
scale at which the Standard Model gauge couplings become unified. Perhaps 
the main message is that the fiatness of the potential has surprisingly little 
effect on the string cosmology. 

We also discuss global defects, and in particular focus on the case of global 
textures and monopoles. We arrive to a rather surprising conclusion that 
global texture and monopole dynamics is not affected at all by the flatness 
of the potential. This has a simple explanation: the potential is simply not 
flat enough! The 'flatness scale' mo ~ lO'^GeV is large in comparison to the 
scale relevant for cosmology: T ~ leV, so the flat potential suffices to confine 
late textures to the vacuum manifold. This leaves both textures and global 
monopoles as a viable candidate for structure formation with the scales of 
symmetry breaking given by M ~ 1.4 x lO^^GeV and M ~ 1 x lO^^GeV, 
respectively. Finally in passing we give an argument that domain walls are 
still a problem. 

In summary, the study of fields with fiat potentials and large vevs is be- 
coming another field on which particle physics increasingly meets cosmology. 
Our hope is that cither cosmology will put constraints on particle physics 
models, or particle physics will offer some interesting cosmological phenom- 
ena, at best tell us something about the origin of structures in the Universe. 
Since this area of research is developing very fast, it is very likely that in near 
future one will be able to make more definite statements about the formation 
and nature of defects in theories with fiat potentials, and hence make a more 
definite prediction on their cosmological relevance. 
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